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Abstract. We present the first high resolution X-ray spectrum of the bright 04Ief supergiant star £ Puppis, 
obtained with the Reflection Grating Spectrometer on-board XMM-Newton. The spectrum exhibits bright emission 
lines of hydrogen-like and helium-like ions of nitrogen, oxygen, neon, magnesium, and silicon, as well as neon-like 
ions of iron. The lines are all significantly resolved, with characteristic velocity widths of order 1000 — 1500 km s~ 1 . 
The nitrogen lines are especially strong, and indicate that the shocked gas in the wind is mixed with CNO-burned 
material, as has been previously inferred for the atmosphere of this star from ultraviolet spectra. We find that 
the forbidden to intercombination line ratios within the helium-like triplets are anomalously low for N VI, O VII, 
and Ne IX. While this is sometimes indicative of high electron density, we show that in this case, it is instead 
caused by the intense ultraviolet radiation field of the star. We use this interpretation to derive constraints on 
the location of the X-ray emitting shocks within the wind that are consistent with current theoretical models for 
this system. 
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1. Introduction 

The initial discovery of X-ray emission from O stars with 
the Einstein Observatory in the late 1970s (Harnden et al. 



197S ) , sparked a vigorous field of research aimed at bet- 



ter understanding the production of hot gas in such sys- 
tems. Most current models invoke hydrodynamic shocks 
resulting from intrinsic instabilities in the massive, radia- 
tively driven winds from these stars (for a detailed review 
see Feldmeier et al. 1997). Support for this picture comes 



from the fact that the X-ray flux is not highly absorbed 
at low energies (Cassinelli & Swank 1983; Corcoran et al. 



1993), which is expected if the X-ray emitting gas is dis- 



tributed throughout the wind rather than in some form of 
hot corona in the outer atmosphere of the star. However, 
the X-ray observations to date have not been especially 
constraining for stellar wind models. This is primarily due 
to the low spectral resolution of the available nondisper- 
sive detectors, which has precluded the study of individual 
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atomic features so crucial to the unambiguous determina- 
tion of physical conditions in the shocked gas. 

In this Letter, we present one of the first high res- 
olution X-ray spectra of an early-type star, the 04Ief 
supergiant £ Puppis, which we have obtained with the 
Reflection Grating Spectrometer (RGS) experiment on the 
XMM-Newton Observatory. £ Pup is an excellent target 
for such work since it is the brightest O star in the sky, and 
has consequently been very well studied at longer wave- 
lengths. In particular, Pauldrach et al. (1994) have devel- 
oped a very detailed NLTE model of the atmosphere and 
wind of this star, constrained by high quality ultraviolet 
spectra from Copernicus and IUE. They find a mass loss 
rate of ~ 5.1 x 10" 6 M© yr _1 and an effective temperature 
~ 42 000 K. The derived abundances indicate that the at- 
mosphere is mixed with CNO-processed material, consis- 
tent with the theoretical picture of a highly evolve d star 
near the end of core hydrogen burning. Hillier et al. ( 1993 ) 
found that the Pauldrach et al. model is consistent with 
the ROSAT X-ray spectrum of ( Pup, but only if the X-ray 
emission arises in shocks distributed throughout the wind. 
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At the very lowest energies (< 200 eV), they find that the 
emitting material must be > 100R* away from the star. 
This is a consequence of the fact that in the Pauldrach 
et al. model, helium is mostly only singly ionized in the 
outer regions of the wind, so that the photoelectric opacity 
is very high at soft X-ray energies. 

Our XMM-Newton RGS spectrum is dominated by 
broad emission lines of mostly hydrogen-like and helium- 
like charge states of nitrogen, oxygen, neon, magnesium, 
and silicon, and neon-like ions of iron. The data provide 
a number of important constraints on the nature and lo- 
cation of the X-ray emitting material in the £ Pup wind. 
Of particular interest is the fact that we see a suppression 
of the forbidden line and enhancement of the intercombi- 
nation line in the helium-like triplets of nitrogen, oxygen, 
and neon. We show that this is a natural consequence 
of the intense ultraviolet radiation field in the wind, and 
that it allows us to place constraints on the location of 
the X-ray emitting shocks relative to the star. Additional 
constraints come from the emission line velocity profiles. 
We show that the data are consistent with the predictions 
of the Pauldrach et al. model. 

In section 2, we describe the details of the RGS obser- 
vation, the nature of our data reduction and analysis, and 
the key observational features of the spectrum. In section 
3, we consider the implications of our results for our un- 
derstanding of £ Pup, and of the X-ray emission from O 
stars in general. 

2. Observations and Data Analysis 

The RGS covers the wavelength range of 5 to 35 A with a 
resolution of 0.05 A, and a peak effective area of about 
140 cm 2 at 15 A. £ Pup was observed for 57.4 ks on 
2000 June 8. The data were processed with the XMM- 
Newton Science Analysis Software (SAS). Filters were ap- 
plied in dispersion channel versus CCD pulse height space 
to separate the spectral orders, and the source region was 
separated with a 1' spatial filter. The background spec- 
trum was obtained by taking events from a region spa- 
tially offset from the source. The wavelengths assigned to 
the dispersion channels are based on the pointing and ge- 
ometry of the telesco pe an d are accurate to ~ 0.008 A 
(see den Herder et al. 2000; this vol.). The effective area 
was simulated with a Monte Carlo technique using the re- 
sponse matrix of the instrument and the exposure maps 
produced by the SAS. Based on ground calibration, we 
expect the uncertainty in the effective area to be less than 
10% above 9 A and at most 20% for shorter wavelengths 



(den Herder et al. 2000). A fluxed spectrum, corrected for 
effective area, with the two first order spectra added to- 
gether to maximize statistics, is presented in Fig. There 
is a small discontinuity in the spectrum near the nitrogen 
Lyman (3 line at 20.91 A caused by a gap between two 
CCD chips. The spectrum was also extracted for analysis 
in XSPEC using standard SAS routines. 

Due to a large solar flare event that occurred close to 
the time of the ( Pup observation, all three EPIC detectors 



and the Optical Monitor on XMM-Newton were switched 
off while this source was observed. Thus, only data from 
RGS are available for analysis. 

2.1. Light Curve 

Previous ROSAT obervations of £ Pup (Berghofer et 



al. 1996 ) show evidence for a 16.67 h period of variability 
in the X-ray emission between 0.9 and 2.0 keV. We have 
extracted a light curve (Fig. |J) of the events in this energy 
range in the RGS data set. There is no significant observed 
variation on this time scale apparent to the eye, and the 
fit to a constant intensity yields xl = 0-9- However, our 
upper limit to the percentage variation at this period is 
not inconsistent with the Berghofer et al. detection. 

2.2. Emission Line Intensities and Emission Measure 
Analysis 

The spectrum displayed in Fig. [I] is composed almost en- 
tirely of emission lines, with a weak underlying continuum. 
The lines are all broadened, with characteristic velocity 
widths of order 1000 km s . In nitrogen, oxygen, and neon 
the He-like forbidden lines are very weak and the intercom- 
bination lines are bright. The line fluxes (Table P were 
measured by taking the integrated flux of the spectrum 
over the line and subtracting a corresponding amount of 
continuum flux. Lines that originate from the same ion and 
that are blended were evaluated as one complex (for exam- 
ple, the He-like triplets or the Fe XVII emission around 15 
or 17 A). The continuum strength was determined by tak- 
ing the flux of a spectral region free of lines but near the 
line in question. In cases where this was not possible, the 
continuum strength was interpolated from the strength 
in other regions of the spectrum. Some of the measure- 
ments were complicated by the presence of overlapping 
lines. When possible, the flux of these other lines was es- 
timated from the flux of lines originating from the same 
ion by comparing with line power ratios. 

We calculated the emission measure (see Fig. ||) for 
each ion assuming solar abundances (Anders & Grevesse 



1989) and a temperature given by the temperature of for- 



mation for the dominant lines from that charge state. The 
emission measure is 

EM = -Fimc Tr~~~~7T (1) 

where -Fii no is the observed flux in the line, d is the dis- 
tance to the source, Pii nc is the line power, A is the ele- 
mental abundance, and /i is the ion fraction evaluated at 
the temperature of formation. We used line powers from 
the APEC code (Smith & Brickhouse |2000|) which in- 
cludes ion fractions from Mazzotta et al. ( 1998| ). We take 



d = 450 pc (Schaerer et al. |1997| ). The line fluxes must 



be corrected for interstellar absorption, although this is 
a min or eff ect. We take iVn = 10 20 cm~ 2 (Chlebowski 
et al. 1989) and we use cross sections from Morrison & 



McCammon (1983) 
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Fig. 1. First order background subtracted spectrum. It has been corrected for the effective area of the instrument. 



We find that the emission measures derived from the 
nitrogen emission lines are at least an order of magnitude 
greater than those for carbon and oxygen, which both 
have temperature of formation ranges that overlap with 
that of nitrogen. This indicates that the ratios A-^/Ao 
and An /Ac are substantially higher than solar, even al- 
lowing for a factor of two uncertainty due to the crude- 
ness of the emission measure analysis. This result is con- 
sistent with the inference of atmospheric abundances by 



Pauldrach et al. (1994), based on their analysis of UV 
spectra. It indicates that the wind material in £ Pup has 
been significantly mixed with matter that has undergone 
CNO burning in the stellar interior. 

The emission measure expected for a smooth, spher- 
ically symmetric wind with parameters appropriate to 
C Pup is EM = 6.5 x 10 60 cm~ 3 . The fact that we find 
EMs that are lower by 4 to 5 orders of magnitude implies 
that only a small fraction of the wind material is heated 
to X-ray emitting temperatures. 
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Fig. 2. Light curve of first order events in the range 6.2 
to 13.8 A, corresponding to 0.9 to 2.0 keV. Events from 
both instruments have been binned together. The time 
dependence is plotted as a function of the 16.67 hr period 
reported in Berghofer et al. 



2.3. Continuum emission analysis 

The continuum emission has a very low intensity relative 
to the bright line emission. To ensure that the contin- 
uum that we see in the spectrum is real, and not, for 
example, an artifact of faulty background subtraction, we 
plotted the first order events in a region free of emission 
lines as a function of their spatial distribution in the cross- 
dispersion direction, without subtracting the background. 
The peak at the location of the source spectrum is clearly 
visible, and indicates the presence of true continuum, well 
above background. The detected continuum also signifi- 
cantly exceeds the scattered light contribution from the 
ensemble of detected emission lines. 



We looked for evidence of discrete photoelectric ab- 
sorption edges in the continuum. Although there are no 
edges evident, the constraints are weak due to the low 
strength of the continuum emission. The upper limits on 
the optical depths are given in Table § These upper lim- 
its are incompatible with the detection of a strong edge 
feature near 0.6 keV reported by Corcoran et al. (1993), 
but that is perhaps not surprising, considering the com- 
plexity of the spectrum and the low spectral resolution of 
their measurements. 

The continuum is so weak that it is difficult to ascer- 
tain its overall shape, especially in the region from 8 - 
17 A where there are many lines. We tried fitting with 
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Fig. 3. Emission measure calculated from each line as a 
function of Tf or m- The horizontal error bars indicate the 
temperature range over which Pine/i is at least half its 
maximum value. Solar abundances are assumed. 

Table 1. Measured fluxes for prominent emission lines 
with rest wavelengths in A. 



Line 


Flux 






(photons cm -2 s~ 


-4 



Mg XII (8.42) 
Mg XI (9.17,9.23,9.31) 
Ne X (12.13) 
Ne IX (13.45, 13.55, 13.70) 
Fe XVII (15.01, 15.26) 
Fe XVII (16.78, 17.05, 17.10) 
O VIII (18.97) 
O VII (21.60,21.80,22.10) 

N VII (24.78) 
N VI (28.78, 29.08, 29.53) 
C VI (33.74) 



(4.1 ± 2.8) x 10~ 5 
(2.0 ±0.7) x 10~ 4 
(1.8 ±0.7) x 10~ 4 
(5.0 ± 1.0) x 10~ 4 
(7.4 ±0.7) x 10~ 4 
(4.7 ±0.8) x 10" 4 

(3.50 ±0.22) x 10" 4 
(5.1 ±0.4) x 10" 4 
(6.4 ±0.5) x 10" 4 

(1.13 ±0.11) x 10~ 3 
(6.4 ±2.8) x 10~ 5 



Table 2. Upper limits on the strengths of the K-edges of 
Ne, O, and N. 



Edge 


Wavelength 

(A) 


Upper limit on To 


Ne X 


9.10 


0.5 


Ne IX 


10.37 


0.5 


O VIII 


14.23 


0.2 


O VII 


16.78 


0.3 


N VII 


18.59 


0.5 


N VI 


22.46 


0.5 



a bremsstrahlung model in XSPEC, but the results were 
inconclusive. 



2.4. He- 1 ike triplet ratios 

The forbidden to intercombination line ratios, R = f /i, 
for helium-like oxygen and neon were obtained by fitting 
in XSPEC. In each case, we fit a gaussian to the Lya 
line of the respective element, and used that profile for 
each of the three components in the He-like triplet. We 
also take into account both the low level continuum and 
lower flux lines in the near vicinity. For Ne IX we find R = 
0.34 ±0.11 and for O VII we find R = 0.19 ± 0.08. The fit 
to Ne IX is shown in Fig. ||. These values are well below 
the expected values for low density plasmas in collisional 
equilibrium (see Mewe et al. 2000). A similar analysis did 
not work very well for the N VI He-like triplet because 
of the complexity of both the Lya and He-like profiles 
for that element. Nevertheless, it is clear from the data 
that the forbidden line is strongly suppressed for N VI as 
well. The Si XIII and Mg XI triplets are too blended to 
allow the forbidden line and intercombination line to be 
quantitatively separated. Again, it is clear from the data 
that the Mg XI forbidden line is suppressed, although not 
as strongly as for Ne IX. 

The conversion of forbidden line to intercombination 
line emission in high density plasmas is a well known ef- 
fect. However, this can also occur at much lower densi- 
ties if the plasma is exposed to a strong UV radiation 
field (Mewe & Schrijvcr 1978). To produce a low R ra- 
tio, electrons populating the 2 3 Si state must be excited 
into the 2 3 P state. The state 2 3 Si is metastable, but 
it is easily populated at collisional equilibrium temper- 
atures, and the corresponding emission line intensity is 
normally strong. When the excitation rate to the 2 3 P lev- 
els becomes comparable to the decay rate to the ground 
state, the forbidden line is suppressed. In a high density 
plasma, this occurs because the collisional excitation rate 
from 2 3 Si to 2 3 P competes effectively with radiative de- 
cay. However, since this is a dipole transition, it can also 
be photoexcited if the ambient UV flux at the appropriate 
wavelength is sufficiently high. To calculate the photoex- 
citation rate, we estimated the emission from £ Pup for 
the frequencies of the 2 3 Si to 2 3 P transitions for N, 
O, Ne, Mg, and Si. We assumed a blac kbod y spectrum 
with T eS = 42 000K (Pauldrach et al. |l994| ). The rate 

of photoionization is given by Ppe = P^t^t/j where / is 
the oscillator strength. We used oscillator strengths from 



Cann and Thakkar ( |1992| ) and Sanders and Knight fll989p . 
The flux is given by F u = 2n(l - yT - (^-) 2 )I V , where 
I v is the specific intensity for the blackbody. 

We used the archived IUE UV spectrum of £ Pup in 
addition to the Copernicus UV spectrum (Morton and 
Underhill 1977) to assess the validity of our blackbody 
model for the UV emission from the photosphere. We com- 
pared the measured flux at the wavelengths of the 2 3 S\ 
to 2 3 P transitions for N VI, O VII, and Ne IX to the flux 
predicted by the blackbody model. After correcting the 
measured flux for absorption, these agree to within a fac- 
tor of two. Since the critical radius where Ppe = Pdecay 
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Table 3. Comparison of photoexcitation and decay rates 
of the 2 3 S state. Decay rates are from Drake (1971). 



13.6 13.7 
wavelength (ii) 



Fig. 4. The Ne IX triplet with best fit model. Only RGS 
2 data are shown for clarity, but both were used in the 
fit. There is also a Fe XVII line at 13.825 A. The rest 
wavelengths of the three lines are at 13.45 A (r), 13.55 A 
(i), and 13.70 A (f). 



depends approximately on \fF\, : the radii we calculate are 
valid to within a factor of The measured UV spectra 
also indicate that there is negligible optical depth in these 
lines due to the wind. This is expected since these are ex- 
cited state transitions, and since the helium-like plasma 
only represents a very small fraction of the wind. 

In Table || we list the decay rates of 2 3 Si and the 
photoexcitation rates from 2 3 Si to 2 3 P evaluated at the 
photosphere. It is clear from this calculation that the for- 
bidden line suppression we observe is in fact due to pho- 
toexcitation from £ Pup's high UV flux, and not due to 
collisional excitation at high densities, as long as the emit- 
ting regions are close enough that the UV is not sufficiently 
diluted. We can thus place constraints on the location of 
shock formation for each of the observed lines. Since the 
forbidden lines in N, O, and Ne are strongly suppressed, 
the emission must occur at radii smaller than the radii at 
which the photoexcitation and decay rates are equal (the 
critical radius). Since Mg XI has a somewhat suppressed 
forbidden line, its emission probably occurs near the crit- 
ical radius. The emission from Si must occur farther out 
than the critical radius. 



2.5. Line profile analysis 

The projected velocity (v p ) profile for a thin spherical shell 
(with a single radial velocity) is flat, or J^- = (const.). 
We expect the emission line profile to appear as a con- 
volution of the radial emission intensity with this flat 
projected velocity profile. The velocity of each shell is 
given as a function of radius by the conventio nal 8- 
model: v(r) ss Wqo(1 — r /r)P (Lamers & Cassinelli 1999 ) 
where /3 « 0.8 and ro ~ i?*. Lines formed at larger radii 
will therefore appear broader than lines formed close to 
the star. Furthermore, lines originating from larger radii 



Ion 


Decay Rate 


Photoexcitation 


Radius 






Rate (R = R«) 


at which 




2 3 S -> 1 \S 


2 3 S -» 2 3 P 


RpE = -Rdecay 








(R») 


Si XIII 


3.56 x 10 s 


4.83 x 10° 


2.7 


Mg XI 


7.24 x 10 4 


7.36 x 10 6 


7.1 


Ne IX 


1.09 x 10 4 


1.11 x 10 7 


22.6 


VII 


1.04 x 10 3 


1.66 x 10 r 


89 


N VI 


2.53 x 10 2 


1.99 x 10 7 


198 



Table 4. Velocity widths and shifts of the Lyman a lines. 



Ion 


Velocity shift 


Velocity width 




(kms -1 ) 


(kms" 1 ) 


Ne IX 


250 ± 125 


940 ± 150 


O VII 


400 ± 80 


1230 ± 80 


N VI 


0±60 


1370 ± 100 



are formed over a region with a small velocity gradient. 
Therfore, we expect these lines to appear more flat-topped 
than lines originating closer to the star. This is apparent 
in the observed Lya lines, which are plotted in velocity 
space in Fig. ||. The N VII peak is noticeably broader and 
has a substantially different shape than the other peaks 
in the plot. It also shows evidence for resolved, discrete 
structure, given the resolution of the instrument. Note 
that other lines overlap with the Ne X line; it does not 
have a red shoulder. 

In Table | we list the shifts in the line centroids and 
the line widths in velocity space for Ne X, O VIII, and N 
VII. The fit to N VII is poor due to the complex structure 
in the line profile. 

3. Discussion 

The RGS spectrum we have presented provides very 
strong confirmation for a number of key aspects of conven- 
tional models of X-ray emission from early-type stars in 
general, and for £ Pup in particular. The emission is dom- 
inated by broad emission lines indicative of hot plasma 
that is flowing outward with the wind. We do not see ev- 
idence for strong attenuation at low energies, which con- 
firms that the X-ray emitting regions are not concentrated 
at the base of the wind, but instead are distributed out to 
large radii. These very simple results agree well with the 
predictions of wind instability models. 

As we have shown, the suppression of the forbidden 
lines in the He-like triplets of low Z elements allows us to 
derive upper limits to the radii of the emitting shocks in 
each case. It is interesting to compare these values with 
Fig. 1 of Hillier et al. ( 1993 ), which is a plot of the radius 
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homogeneity. With its large projected rotational velocity 
of v sin i ~ 203kms _1 (Penny 1996), CP U P most probably 
falls into this category. 
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Fig. 5. Ly a lines in velocity space. Nitrogen is the solid 
line, oxygen is dashed, and neon is dotted. The intensities 
have been renormalized for comparative purposes. The er- 
ror bar is representative for the nitrogen line. Note the 
discrete structure in the nitrogen peak. The shoulder on 
the Ne X line is not a velocity feature; it is an emission 
line from Fc XVII. 



at which optical depth unity is reached in the wind as a 
function of X-ray energy for th e NL TE model of £ Pup 
calculated by Pauldrach et al. ( 1994 ). At the energies of 
the N VI, O VII, Ne IX, Mg XI, and Si XIII lines, unit 
optical depth is achieved at 22, 23, 9, 3.5, and 2.5 stellar 
radii, respectively, in this model. These values are quite 
compatible with both our derived upper limits for N VI, 
O VII, Ne IX, and Mg XI, and our derived lower limit for 
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